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ABSTRACT 

Context. The A-variance analysis is an efficient tool for measuring the structural scaling behaviour of interstellar turbulence in 
astronomical maps. It has been applied both to simulations of interstellar turbulence and to observed molecular cloud maps. In paper 
I we proposed essential improvements to the A-variance analysis and tested them on artificial structures with known characteristics. 
Aims. In this paper we apply the improved A-variance analysis to simulations of interstellar turbulence and observations of molecular 
clouds. We tested the new capabilities in practical use and studied properties of interstellar turbulence that could not have been 
addressed before. 

Methods. We selected three example data sets that profit in particular from the improved A-variance method: i) a hydrodynamic 
turbulence simulation with prominent density and velocity structures, ii) an observed intensity map of p Oph with irregular boundaries 
and variable uncertainties of the different data points, and iii) a map of the turbulent velocity structure in the Polaris Flare affected by 
the intensity dependence on the centroid velocity determination. 

Results. The tests confirm the extended capabilities of the improved A-variance analysis. Prominent spatial scales were accurately 
identified and artifacts from a variable reliability of the data were removed. The analysis of the hydrodynamic simulations showed 
that the injection of a turbulent velocity structure creates the most prominent density structures are produced on a scale somewhat 
below the injection scale. The new analysis of ap Oph continuum map reveals an intermediate stage in the molecular cloud evolution 
showing both signatures of the typical molecular cloud scaling behaviour and the formation of condensed cores. When analysing the 
velocity structure of the Polaris Flare we show that a universal power law connects scales from 0.03 pc to 3 pc. However, a plateau in 
the A-variance spectrum around 5 pc indicates that the visible large-scale velocity gradient is not converted directly into a turbulent 
cascade here. It is obvious that, for any turbulent structure, effects of low-number statistics become important on the driving scale. 

Key words. Methods: data analysis - Methods: statistical - ISM: clouds - ISM: structure 



1. Introduction 

Observations of interstellar clouds show a complex, filamentary 
structure which can be attributed to turbulence in the inte r stellar 
medium (Franco & Carrami nanaL 1 1999b lOssenkopf et all [2000; 
Mac Low & Klessen, 2004). To understand the processes gov- 
erning the structure and evolution of the clouds, turbulence mod- 
els have to be constructed and compared to observational data. 
Their parameters and implementational details need to be ad- 
justed to fit the observed behaviour. Due to the random nature 
of turbulence, simulations will never provide an exact reproduc- 
tion of the observed data sets but will only reproduce general 
statistical properties like scaling relations. 

The most common scaling relation to characterise turbulent 
structures is the power spectrum of fluctuations, both in den- 
sity and in velocity. Here, we encourage using another quantity, 
the A-variance spectrum, introduced by Stutzki et al.l dl99 8). a 
wavelet-based method to measure the relative amount of struc- 
tural variation as a function of the size scale. Due to the mutual 
relations between the A-variance spectrum and the radially aver- 
aged power-spectrum, the A-variance analysis can be considered 
as a very robust method of evaluating the power spectrum of a 
structure. The advantages of the A-variance method result from 
the smooth wavelet filter shape, which provides a robust way for 



an angular average independent of gridding effects, and from th e 
insensitivity to edge effects as discussed by Bens ch et al.l (12001). 



In parallel to the structure scaling analy sis, clump decom- 
positi on algorithms like GAUSSCLUMP dStutzki & Giistenl 
1990) found that the clump mass spectrum dN/dM oc M~ y 
of sufficiently large molecular cloud data sets also tend to 
follow power laws over many orders of magnitude with the 
spectral index y in a relatively narrow ra nge betw een 1.7 to 
1.9 dKramer et all 119981: iHeithausen et all 1 1998b . IStutzki et at] 
(1998) demonstrated that a clump ensemble with such a mass 
spectrum and with a power-law mass-size scaling relation re- 
sults in a cloud image with the A-variance scaling index deter- 
mined by the number-mass and the mass-size spectral indices. 
However, it was questioned by e.g. [V azquez-Semaden i et all 
dl997l) ; lBallesteros-Paredes & Mac Low! d2002l) whether the ob- 
served mass-size relation reflects true properties of the underly- 
ing structure. It may rather represent an observational artifact. 
From these results it is obvious that there is a strong interest in 
studying the spatial scaling behaviour of observed maps of the 
interstellar medium via the A-variance analysis. 

In paper I we have proposed several essential improvements 
to the original A-variance method. We have investigated the 
use of different wavelets and calibrated their spatial resolution. 
Unfortunately, it turns out that it is not possible to define a sin- 
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gle optimum wavelet for all purposes because different wavelets 
exhibit a different power in the detection of the characteristic 
structures. A good compromise is given by the Mexican-hat fil- 
ter with a diameter ratio of 1.5. When the main focus lies on the 
measurement of the spectral index, the French-hat filter with a 
diameter ratio of about 2.3 is also suitable. We also introduced 
a significance function to weight the different data points. This 
allows us to analyse observed data where the signal-to-noise ra- 
tio is not uniform across the mapped area, but spatially varying. 
This should permit to distinguish the influence of variable noise 
from actual small-scale structure in the maps. The need for such 
a trea tment became very obvious when Ossenko pf & Mac Low! 
(2002) used the A-variance analysis to characterise the velocity 
structure detected in molecular line obser vations of the Polari s 
Flar e taken by Fal garone et al. (1998), iBensch et ail (1200 ll) , 
and Heithausen & Thaddeus (1990). Comparing the A-variance 
analysis with the size-linewidth relation, we found that the A- 
variance of the centroid velocity maps produced wrong results 
on scales where the maps do not show any noticeable emission. 
With the introduction of the significance function, the new A- 
variance analysis should be able to also reliably analyse such 
data sets. Moreover, the use a weighting function also allows 
us to use computational methods like the fast Fourier transform 
to obtain the A-variance spectrum even for maps with irregular 
boundaries or maps which are only sparsely filled by significant 
values. 

In paper I we tested the properties of the improved A- 
variance analysis using simple artificial data sets. In this paper 
we will apply it to more realistic data sets, either simulations 
of interstellar turbulence or to actual observational data. In Sect. 
2 we recapitulate the formalism of the new A-variance analy- 
sis and summarise the results that we obtained from the appli- 
cation to the test data sets. In Sect. 3 we apply the analysis to 
a hydrodynamic sim ulation, to the p-Oph dust continuum map 
by M otte et all dl99 8) and to the centroid velocity map for the 
Polaris Flare. We discuss the conclusions on their structure in 
Sect. 4. 



2. The improved A-variance method 

In this section we summarise the main properties of the im- 
proved A-variance analysis proposed in paper I, focusing on the 
new points, not cov ered in the original A-variance definition by 
IStutzkiet all (119981) . 

The A-variance measures the amount of structure on a given 
scale / in a map f(r) by filtering the map with a spherically sym- 
metric wavelet of size / and computing the variance of the thus 
filtered map: 
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where, the average is taken over the area of the map, the symbol 
* stands for a convolution, and Q, describes the filter function 
composed of positive inner "core" and a negative annulus, both 
normalised to integral values of unity 
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We have studied the "French-hat" filter with constant values in 
both parts: 
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and a "Mexican hat" consisting of two Gaussian functions: 
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where I is the core diameter and v is the diameter ratio between 
the annulus and the core of the filter. Plotting the A-variance as a 
function of the filter size / then provides a spectrum showing the 
relative amount of structure in a given map as a function of the 
structure size. 

The effective filter size, given by the average distance of 
points in the core and the annulus, deviates from the core di- 
ameter / as 

l eB _ j 0.29v + 0.26 for the French hat 

I " 1 0.41v + 0.46 for the Mexican hat. 



(5) 



Thus structures with a particular size should show up as promi- 
nent peaks in the A-variance spectrum on a scale l^ correspond- 
ing to that size. Test with artificial data sets in paper I have 
shown, however, that the peak positions always falls 10-20% 
below the maximum structure size. Taking this systematic off- 
set into account we can nevertheless reliably calibrate the spatial 
resolution of the A-variance analysis. 

A major improvement of the new A-variance algorithm was 
the introduction of a weighting function to the data w^at/)- This 
simultaneously solved the problems of the edge treatment of fi- 
nite maps and the analysis of data with a variable uncertainty 
across the map. The weight function varies between and 1, 
representing the reliability of the individual data points, and it 
extends beyond the original map size, padding it with zeros at 
the boundaries. Instead of the original map f(r), an extended 
map, /paddedW = f( r ) x Wdata(>) inside the original data area, 
/padded(f) = outside, is analysed. This padded map can be peri- 
odically continued without wrap-around effects, so that the filter 
convolution can be efficiently computed in Fourier space involv- 
ing a fast Fourier transform and a map multiplication. 

To avoid that data points within the padded area or with a 
low weighting are counted like normal zero-value data, but dis- 
regarded in the computation of the variance, the filter has to be 
re-normalised at each position in the map in such a way that the 
integral weights of core and annulus remain unity when exclud- 
ing the padded points and when taking the weighting of the nor- 
mal data points into account. Instead of one convolution (Eq.[TJ, 
one has to compute four convolutions 



G/.coieO) = /paddedO) * QO"') 

/,core 

Gu nn (r) = /paddedW * 0( r ') 
/,ann 

W,, core (r) = w(r) * Q(r') 

/,core 

W) = *f(r) * Q(r') 
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and combine the results while re-normalising with the effective 
filter weight for the valid data 

G/.coref/) Gi tWm (r) 
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From the actual filter weight computed for each point in the 
map we can derive a significance function as the product of both 
normalisation factors 



W,, tot (r) = W,, core (r)W / , ann (r) . 



(8) 



This provides the actual significance of the data points in the 
convolved map which is used when computing the A-variance of 
the whole map 
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With this generalised concept, the A-variance analysis can be 
applied to arbitrary data sets. They must be projected onto some 
regular grid but they do not need to contain regular boundaries as 
the corresponding "empty" grid points can be marked with a zero 
significance. Varying noise or other changes in the data reliabil- 
ity can be expressed in the significance function WdataO")- This 
applies e.g. to maps where not all points are observed with the 
same integration time so that they show a different noise level. In 
paper I we used a weighting function given by the inverse noise 
RMS and in Sect. 13. 21 we will study the impact of the selection 
of the weighting function for observed data. 

The only remaining requirement for the applicability of the 
A-variance analysis is a sufficiently large spatial dynamic range 
in the data. B ensch et alj d200ll) had shown that a map has to 
contain at least 30 pixels in each direction to obtain reasonable 
error bars of the A-variance spectrum. Numerical tests with noisy 
data in paper I showed that this critical size needs to be extended 
by a factor of about one over the average data significance in the 
case of data with a variable reliability. 

3. Applications 

3. 1 . Hydrodynamic simulations 



In the 



papers 



lOssenkopf & Mac 



> by JMac Low & O ssenkopf (|2000h . 

Low! J2002I) . and|Ossenkopfetal. (2001 we 



have demonstrated the general applicability of the A-variance 
analysis to extract characteristic structure sizes and scaling 
laws from (magneto-)hydrodynamic simulations performed 
with a variety of codes. Here, we need to test whether the 
A-variance with adapted filter functions improves the sensitivity 
of this method. The weighting function is irrelevant in this case 
because the data do not suffer from noise or another cause of 
variable reliability across the data set. 

We h ave applied the analysis to a variety of simulations pre- 
sented in lOssenkopf et alj d200ll) . but we present the results here 
only for a single model, the first inertial stage of the small-scale 
driven hydrodynamic turbulence computed by smooth-particle 
hydrodynamics (SPH), S02 at t = 0. In this simulation the ve- 
locity field is driven by a Gaussian field of random fluctuations 
within a finite wavenumber range, k - 7 ... 8. This means, that 
the driving process introduces characteristic variations into the 
velocity structure with the same scale length as used in the ar- 
tificial sine wave field used in paper I, but with wavenumbers 
between 7 and 8. Thus the A-variance spectrum should mea- 
sure a peak variation for scales of l/( V2A:), i.e. between 0.088 
and 0.101 of the size of the whole data cube. Selecting a model 
which is driven on small scales guarantees that we can identify 
a clear peak for these structures leaving enough dynamic range 
on smaller and larger scales. We select the initial stage of fully 
evolved turbulence in the simulation to make sure that the tur- 
bulent driving is the only process creating structures in the data 
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Fig. 1. A-variance spectra of the z-projection of the z-component 
of the velocity field in the hydrodynamic simulation S02 at t = 0. 
Four different filter shapes are tested against each other. 



set, avoiding effects of self-gravity. In this way we have a data 
set which is best suited to test the structure recognition by the 
A-variance analysis, as it should directly detect the scales of the 
driving process in the velocity structure and determine the scal- 
ing of the turbulent cascade on lower scales. Results for the other 
simulations did not provide any fundamentally different results, 
but are less clear to interpret because either the driving scale is 
closer to the edge of the dynamic range or the contained struc- 
ture is less well known. 

Figure [1] shows the A-variance spectra of the z-projection of 
the z component of the velocity field of the turbulence simulation 
computed with four different filter functions. We find a clear shift 
between the peak positions measured by the different filters. The 
Mexican-hat filter gives systematically larger lags for the peak 
than the French hat and the lag of the peak grows with growing 
annulus-to-core diameter ratio, v. The minimum peak lag, given 
by the French hat with v = 1 .5, falls at 0.084, the maximum lag, 
given by the Mexican hat with v = 3, at 0.11. 

This behaviour is consistent with the results obtained for the 
simple sine wave field in paper I. It can be understood as a re- 
sult of the variable width and the shape of the filters in Fourier 
space. The broadening of the filter function reduces the contrast 
which leads in a spectrum with a steep slope on small scales 
and a shallow slope on large scales to an effective shift of the 
peak position. The steep decay of the French-hat filter function 
for large lags leads to a somewhat lower peak position, but is 
always accompanied by side lobes of the Bessel function visi- 
ble as artificial secondary peaks at large lags. For the turbulence 
simulations these secondary peaks are not as pronounced as for 
the sine wave field, but they are also visible in Fig. FJ For the 
diameter ratios v of about 1.5 for the Mexican hat and 2.3 for 
the French hat, deduced as optimum values in paper I, the peak 
position falls at about 0.088 in both cases, a value slightly lower 
than the expected average structure size. 

The different width of the peak has a direct impact on the 
slope of the turbulent structures measured at small lags. With 
the broad peaks produced by the Mexican-hat filter, the slope 
is affected down to relatively small scales. As the filter diame- 
ter ratio v constrains the minimum scale which can be resolved, 
a clear power law becomes only visible for the French-hat fil- 
ter with v = 1.5. Small diameter ratios are always favourable 
with respect to the dynamic range which can be covered in the 
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Fig. 2. A-variance spectra of the density structure created by 
the velocity field from Figfl] measured with four different filter 
functions. 



A-variance analysis because of the minimum filter size and the 
constraint that the overall filter must always remain small com- 
pared to the analysed map. For the French hat with v = 3.0 and 
the Mexican hat with v = 1 .5 the remaining dynamic range is 
just marginally sufficient to reliably determine the spectral in- 
dex. We find slopes between 2.8 and 3.1 corresponding to power 
spectral indices £ - 4.8 ...5.1. This is much higher than the 
Kolmogorov index of £ = 3.67 indicating that the scaling on 
small scales is not determined by a self-similar turbulent cas- 
cade, but by the numerical viscosity in these simulations damp- 
ing small scale structures. For simulations driven on large scales 
resulting in a larger dynamic range below the peak, we find a 
limited inertial range with a slope of about two, corresponding to 
£ » 4 which is consistent with a cascade of Burger's t urbulence 
(fossenkopf & Mac LowLl2"002tlOssenkopf et alffcooU see). 



In a second step we investigate how the turbulent velocity 
scaling translates into the creation of turbulent density enhance- 
ments. Figure [2] shows the A-variance spectra of the density 
structure seen in the same step of the simulation using the same 
filters as applied to the velocity structure. We find the same sys- 
tematic deviations between the results seen by different filters 
but a generic shift of the peak position to shorter lags by a factor 
0.75-0.8 with respect to the peaks of the velocity structure. The 
slope of the turbulent density structure on small scales is shal- 
lower by 1 .0 + 0.05. For the mutual comparison between density 
and velocity structure, the selection of the filter is thus irrelevant 
as long as the same filter is applied in both cases. 

We see that injesting energy on a particular scale does not 
create density enhancements on that scale, but rather on a scale 
smaller by a factor 0.75-0.8. Thi s shift has not been noticed be- 
fore by iMac Low & Ossenkopi] (|2000) as only the systematic 
tests of the filter functions provided enough sensitivity with re- 
spect to a reliable scale detection. It seems that the turbulent cas- 
cade builds up density fluctuations on all scales below the driv- 
ing scale, but that those density enhancements act themselves as 
points of an efficient energy conversion between the scales cre- 
ating new density structures so that dominant density scale falls 
somewhat below the initial scale. 

Finally we need to address the significance of the mea- 
sured structure size and the scaling indices with respect to 
the random fluctuations always present in turbulence simu- 
lations and with respect to the relation between the three- 
dimensional (3-D) structure and the two-dimensional (2-D) pro- 



Fig. 3. A-variance spectra of the different rectangular projec- 
tion maps and of the full three-dimensional data cube of the z- 
component of the velocity field in simulation S02 at t = 0. The 
analysis used the French-hat filter with a annulus-to-core diam- 
eter ratio v = 2.3 



jections which we can measure in astronomical observations. For 
turbulent density structures this comparison has been done by 
IMac Low & Ossenkopi] (2000) using the "traditional A-variance 
filter with the diameter ratio v = 3.0. We will repeat it here for 
the velocity structure as the direct carrier of turbulent energy us- 
ing one of more sensitive filter functions. 

Figure [3] shows the A-variance spectra of the z-component of 
the hydrodynamic simulation introduced above, computed with 
a French-hat filter with v = 2.3. The three broken lines show the 
A-variance spectra computed for the three different orthogonal 
projections of the velocity data cube. The solid curve shows the 
spectrum computed for the 3-D structure but rescale as if com- 
puted in 2-D by a factor proportional to the lag, to compensate 
for the different exponents of the A-variance spectra depending 
on the dimensionality of the considered space, and shifted by 
a factor n/4, to account for the average reduction of a random 
structure size when projected from 3-D onto a plane. The solid 
line, thus represents our best knowledge on the velocity struc- 
ture actually present in the simulations, while the broken lines 
represent possible observer's views onto that structure. 

We find a considerable variation between the A-variance 
spectra seen in the different directions, giving a feeling for the 
statistical uncertainty when measuring the scaling in turbulent 
simulation^. Nevertheless, both the peak position and the ex- 
ponent at small lags agree between all three curves. Compared 
to the full 3-D structure, there is, however, a systematic shift of 
the peak to somewhat larger lags and a slight reduction of the 
slope at small lags. For the 3-D structure the peak is seen at a lag 
of about 0.075, which is about 20 % smaller than the expected 
scale for the maximum variation, while the peak for the projec- 
tions falls between 0.081 and 0.092, i.e. only 10 % below the ex- 
pected scale. In the 2-D projections, the peak is always broader 
and the slope at small lags is always somewhat shallower, similar 
to the impact of broader filter functions. It is important to notice, 
that the corresponding plots for the turbulent density structure, 
showing a shallower scaling at small lags, exhibit a very good 
match between the A-variance spectra computed in 2-D and in 



1 The computational uncert ainty given by the fi nite size of any data 
set was discussed in detail by Bens ch et al.l i200lh . 
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3-D, i.e. neither a shift of the peak nor different slopes at small 
lags. 

This seems to indicate that the turbulent velocity cascade 
does not behave fully isotropic. A similar effect would be ex- 
pected for grid-based hydrodynamic simulations where the dis- 
sipation is strong between neighbouring cells in the x-, y-, and 
z-directions. However, the simulations studied here used an SPH 
code which should not show any intrinsic anisotropy. The cause 
for the anisotropy of the turbulent velocity scaling is thus so far 
unknown. 

3.2. Maps with variable noise 

Maps with variable noise are obtained e.g. in observations with 
detector arrays showing a pixel-to-pixel variation in the sensitiv- 
ity. They are produced in single-pixel observations when a drift 
in the receiver sensitivity or the atmospheric conditions changes 
the noise in the data during the measurement and they result from 
mosaicing observations with variable integration times for differ- 
ent regions of the field. All these cases can be analysed in terms 
of the improved A-variance as long as the spatial distribution of 
the noise across the map is well known so that a correspond- 
ing significance map can be defined which is used to weight the 
different points in the A-variance analysis. Then the A-variance 
spectrum is able to distinguish between small-scale noise con- 
tributions in regions with a high noise level and real small-scale 
structures in regions with a low noise level. 

As a challenging example of a data set where a variable noise 
is produced by the observation of different points of the map with 
varying integration time we use the 1 .3 mm continuum map of 
p Oph obtained by Mot te et ail (QJ998). It is the result of a mo- 
saicing observation where an efficient use of the array receiver is 
given by several observations of the source with different orien- 
tations of the array. The combination of these observations then 
results in a poorer coverage of the outer regions of the source 
compared to the central regions in terms of the total integra- 
tion time spent on each point. If the source is covered in total 
with Mot observations, we can characterise the integration time 
at each point by the number of coverages including this point 
1 < iV < iVmax, where Af max < N tot . As the noise at each point is 
inversely proportional to the square root of the integration time, 
we can use the value of V^V as a measure for the data reliability 
across the map. 

In addition to the variable noise the map has highly irregular 
boundaries. Figure [4] shows the intensity map. In contrast to the 
original publication we show the intensity with equidistant con- 
tours on a linear scale because the linear-scale presentation gives 
a better feeling for the structure that is measurable by means of a 
statistical method like the A-variance analysis. To emphasise the 
irregular noise behaviour we plot one contour at 20 mJy/15"- 
beam, which is below the noise level in the outer parts of the 
map and above the noise level in the inner parts. Consequently, 
this contour shows partly real structure and partly artificial struc- 
ture from the noise. Figure [5] contains the corresponding map of 
significance values defined as the square root of the number of 
integrations at each point, yN, thus measuring the inverse noise 
RMS. 

To demonstrate the influence of the significance weighting 
we show in Fig. [6] the A-variance spectra computed for the p 
Oph map using three different weighting functions. The lower 
spectrum is generated when the weights are ignored, i.e. simply 
set to unity at all valid data points. The upper graph is produced 
when the full weighting function from Fig. [5] is used. The in- 
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Fig. 4. 1.3 mm continuum map of p Oph taken by iMotte et al.l 
(199 &). The in tensity is normalised to a 15" beam. In contrast 
to lMotte et al.l we plot the intensity on a linear scale because the 
structure analysis by means of the A-variance only measures this 
linear behaviour. For a better contrast of the picture the intensity 
scale is truncated here at 1000 mJy/beam whereas the map con- 
tains a few points with intensities up to 1700 mJy/beam in the 
brightest core. 




Fig. 5. Map of weights for the intensity map from Fig. |4] given 
by the square root of the number of integrations at each point. 



termediate curve follows when we introduce and upper limit to 
the weighting function motivated by the idea that above a cer- 
tain significance limit a further reduction of the noise level does 
not improve the structure characterisation any more. All com- 
putations use the filter truncation outside of the area where data 
have been taken. The irregular boundaries make it impossible to 
construct a useful periodic continuation here. 

The A-variance spectrum computed without weights seems 
to indicate a wide range of scales with a power-law behaviour 
from about 0.2' to 2' whereas the A-variance spectrum com- 
pute d with the ful l weig hts shows a steepening starting at 0.5- 
0.7'. lBensch et al.l (12001) have shown, however, that such a be- 
haviour is exactly to be expected from the finite beam of the 
observations. The data are given at a resolution of 15" and 
the corresponding beam smearing is known to steepen the A- 
variance spectrum up to scales of about one arcminute. This 
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Fig. 6. A-variance spectra of the 1.3 mm p Oph map computed 
using three different weighting functions. Either all observa- 
tional weights are ignored, they are completely taken into ac- 
count, or weights above y/N = 20 are set to this maximum level. 



steepening can be modelled theoretically using the ana lytic ex- 
pressions for the beam convolution from iBensch et all In fact, 
the A-variance spectrum computed from the p Oph map with 
full weights can be fitted by a single power-law structure with 
a = 0.68 from ^ 0.2' up to about 2' and a convolution func- 
tion of a 15" HPBW beam (see Fig. |7). The fitted exponent of 
0.68 falls into the range measured in molecular line observations 
of molecular clouds covering exponents between 0.5 and 1 . 3 



dBensch et alll200ll:lElmegreen & Scaiotoo^lFafgarone et al.L 
2004). In contrast, the two lower curves cannot be fitted in the 
same way. These spectra would lead to the conclusion of a sur- 
plus of small-scale structure relative to a power-law scaling rela- 
tion. Such a relative surplus of structure on small scales is hard to 
explain as it would require additional driving processes on these 
scales overcompensating the known dissipation of turbulence 
by am bipolar diffusion and molecular viscosity dKlessen et al.L 
2000). Gravitational collapse is not able to create these struc- 
tures; it always affe cts the whole A-varianc e spectrum, not only 
the small-scale tail dOssenkopf et alluOOll) . A surpl us of small- 
scale structure is also in contrast to the analysis of Mott e et all 
dl998l) who found a relative lack of small structures in terms of 
a flatter clump mass spectrum for small clumps. 

Thus we conclude from the scaling behaviour that the full 
weighting of intensity maps by their inverse noise RMS results in 
the most reliable A-variance spectra. With this weighting the A- 
variance analysis is able to distinguish insignificant small-scale 
structure, dominating the lowest contour in Fig. [4] from signifi- 
cant structures which are intuitively be tter presen ted by the con- 
tours chosen in the original plot by Motte et all The increase 
of the absolute value of the A-variance at large lags when us- 
ing the weighting function is explained by the relative increase 
of the contribution of the bright cores in the map when virtu- 
ally reducing the map size by weighting the outer parts by lower 
significance values. 

To get a feeling for the reliability of the different points in the 
A-variance spectrum we plot in Fig. [7] the A-variance spectrum 
including the error bars. The error bars arise from the statisti- 
cal uncertainty of the measurement of the average variance in a 
filtered map (Bensch et al. 2001). Due to the lower number of 
statistically independent points in maps convolved with a larger 
filter, the A-variance is most uncertain at the largest lags. In spite 
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Fig. 7. A-variance spectrum for the p Oph continuum map with 
the full weighting function plotted with error bars indicating the 
statistical uncertainty of the measurement. The connecting line 
represents the power-law fit to the spectrum including the beam 
convolution effect. The point at the smallest lag cannot be fitted 
by the first order approximation of the beam shape (see Bensch 
et al. 2001). 



of the large error bars, the general scaling behaviour can be ac- 
curately traced. The solid line shows the fit to the data using a 
power-law description of the structure scaling and the convolu- 
tion by a 15" beam. 

In the A-variance spectrum one can clearly see that the dom- 
inating structure has a scale of about 2.2', i.e. 0.1 p c. This corre- 
sponds to the typical size of the cores identified by Motte et all 
The contribution of significant structures is continued up to 7-9', 
i.e. 0.3-0.4 pc. This scale agrees with the size of the largest iden- 
tified core. Above about 9' the A-variance spectrum decays with 
a — -2 indicating a lack of further correlated structure on larger 
scales. It is not clear whether the subtraction of large-scale emis- 
sion unavoidable in the used observing mode has removed some 
large-scale correlation which is present in the cloud but cannot 
be detected from the map. 

At small lags the spectrum indicates no separation between 
the scales where "core s" and "clumps" (condensations) were 
defined by iMotte et all The fit in Fig. [7] shows that on scales 
from 0.2' to 2', i.e. within a dynamic range of a factor ten, 
the spectrum is described by a single power-law smoothed by 
the observational beam. The A-variance spectrum suggests that 
the same processes drive the formation of the somewhat larger 
"cores" and the somewhat smaller condensations. The break 
in the spectrum at 0.4', that seems to suggest a change in the 
scaling law of the observed structure, is only produced by the 
beam smearing and is quantitatively in agreement with a con- 
tinuation of the power law observed on larger scales down to at 
least 0.2'. The underlying structure can be described by a per- 
fect a pow er law in co ntradiction to the clump mass spectrum 
studied by IMotte et all who found a significant turn-down at a 
mass of about 0.5 Mq, This difference is even more intrigu- 
ing because of the opposite situation in molecular line studies 
of the Polaris Flar e where the clump mas s spectrum shows a 
perfect power law dHeithausen et all Il998l) but the A-variance 
spectrum shows a steepening towards small scales (Bensch et 
al. 2001, see also Ossenkopf et al. 2000). From the theoret- 
ical modelling of the translation of a clump spect rum into a 
corresponding power spectral index of an fBm by IStutzki et al.l 



V. Ossenkopf, M. Krips and J. Stutzki: A-variance analysis of interstellar turbulence 



( 1998), we would expect a fixed relation between the measured 
clump mass spectrum and the corresponding A-variance spec- 
trum in both cases. However, our examples violate this relation. 
The mass-size relation of the clumps may be affected by opti- 
cal depth effects and a large part of the observationally iden- 
tified clumps may result from the superposition of different 
structures along th e line of sight, not w ell separate d in veloc- 
ity space (Ballesteros-Paredes & Ma c LowL 120021: [Ossenkopf. 
2003). Further systematic studies are necessary to understand the 
actual physical processes interrelating the structure size spectra 
and the clump mass spectra. 

We can compare the new results with the outcome of pre- 
vious A-variance analyses, because all previous conclusions on 
the slopes of the A-variance spectra remain valid. The new A- 
variance method has improved our ability to precisely detect 
prominent scales, it has calibrated the absolute scales and it in- 
creased the statistical significance of the spectra by taking vari- 
able data reliability and edge effects into account, but none of 
these points should significantly affect the general scaling be- 
haviour measured in our previous papers. The p Oph map shows 
a behaviour which i s intermediate between that observed e.g. by 
Bens ch et al.l (|2001 ) in molecular molecular lines, where we find 
a power-law A-variance spectrum on small scales and a dom- 
inance of large-scale structure, and the spectrum measured for 
the 1.3 mm continuum map of Serpens (Testi & Sargent, 1998) 
analysed by Ossenkopf et al. (2001), where small cores domi- 
nate the spectrum resulting in a steep decay on large scales. In 
the p-Oph map we find both effects in one spectrum. The dense 
cores represent the dominating size scale but we can clearly re- 
solve the scaling of significant structure on smaller scales. 

Taking all facts from the scaling behaviour and the clump 
size detection together we conclude that an appropriate charac- 
terisation of the map in terms of a A-variance spectrum is only 
obtained when weighting the intensity maps by the inverse noise 
RMS. Otherwise a variable noise in the spectra always tends to 
mimic small-scale structure which might be taken for real. 

3.3. Velocity centroid maps 

The situation is more difficult in the analysis of maps represent- 
ing other quantities than intensities. Then the weighting function 
given by the inverse noise RMS of the observation is not neces- 
sarily a good measure for the significance of the data. We study 
one such example here. 

Investigating t he velocity structure in the Polaris Flare 
molecular cloud, Ossenkopf & Mac Low (2002) applied the or- 
dinary A-variance analysis to maps of centroid velocities in CO 
line data. The scaling behaviour of the velocity field measured in 
terms of the A-variance spectrum was compared to model sim- 
ulations of interstellar turbulence. The data sets were given by 
three nested CO maps taken with different telescopes at different 
resolutions. The m aps taken at hi gh resolutions with the IRAM 
30 m telescope dFalgarone et all 1 19981) and the KOSMA 3 m 
telescope (Bensch et al. 2001) only covered regions with suf- 
ficiently bright emission so that a reliable determination of the 
line centroid velocities was possible at all points. The centroid 
A-variance spectra derived for these two maps showed a contin- 
uous power-law spectrum with a slight steepening towards the 
smallest lags. 

In contrast, t he map on the largest scale take n with the CfA 
1.2 m telescope (Heithausen & Thaddeus, 1990) contains many 
data points where no emission above the noise limit was de- 
tected. Moreover, it was difficult to obtain a reliable determi- 
nation of the centroid velocities in regions where the line inten- 




Fig. 8. Map of velocity centroids measured in the CO 1-0 
map of the Polaris Flare taken with the CfA 1.4 m telescope 
dHeithausen & T haddeus. 1990). For all points with an inte- 
grated intensity below 0.65 K km/s a reasonable determination 
of the centroid velocity was impossible so that the average ve- 
locity of the cloud (-2.97 km/s) was assigned there. The axes are 
labelled relative to the zero position of / = 123!628, b = 24?93 
in Galactic coordinates. 




Fig. 9. CO 1-0 intensity map corresponding to the velocity cen- 
troid map in Fig. [8] The intensities provide a measure for the 
significance of the centroid velocities. 



sities only exceeded the noise RMS by a factor of a few (see 
Ossenkopf & Mac Low 2002). The resulting A-variance spec- 
trum did not show a continuation of the power-law behaviour 
from the two maps on smaller scales, but turned essentially flat. 
This is in contradiction to an eye inspection of the centroid map 
plotted in Fig. [8] showing a large-scale velocity gradient which 
should appear as well as large-scale structure in the A-variance 
spectrum. The corresponding map of line integrated intensities, 
plotted in Fig. [9] shows that the map contains large regions with- 
out emission. When the ordinary A-variance counts their cen- 
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Fig. 10. A-variance spectra of the Polaris Flare velocity centroids 
computed with different weighting functions characterising the 
significance of the velocities at each point. 



titative impact of the observational noise on the noise in centroid 
velocities. 

In spite of the uncertainty of the noise contribution at small 
lags, we can draw new essential conclusions on the velocity 
structure of the Polaris Flare from the A-variance spectrum. The 
power-law scaling behaviour of the velocity structure detected 
previously only for smaller sizes is now continued up to scales 
of 1 .3°, i.e. 3-4 pc. A single power law with an exponent a « 0.9 
can be used to cover the length scale range of about a factor 100. 
The A-variance spectrum also shows a plateau around 2°, i.e. 
5 pc, indicating the relative deficiency of motions on that scale. 
Above 3° the A-variance spectrum rises again tracing the global 
velocity gradient visible in Fig. [8] This behaviour indicates that 
the large-scale gradient is not converted directly by shear mo- 
tions into the turbulent cascade but that the turbulent cascade 
starts at somewhat smaller scales. This points towards a shock 
producing the large-scale gradient which is only converted into 
turbulent energy on the scale of previously existing density fluc- 
tuations in the cloud and excludes Galactic rotation as the main 
driving force. 



troid velocity with the same weight as that from points in the 
actual molecular cloud the "empty regions" statistically hide the 
variations in the regions with significant values. The virtual lack 
of large-scale velocity variations in the A-variance spectrum is 
thus due to the missing significance weighting. 

The knowledge on the reliability of the centroid veloci- 
ties coming from the corresponding line intensities has to be 
taken into account using the improved A-variance analysis. 
Unfortunately, it is not obvious how the significance of the cen- 
troid velocities is related to the line intensities. We have tested 
weighting functions based on three assumptions: i) The signifi- 
cance of the centroid velocities is determined by the integrated 
line intensities at each point, ii) The zero value of the weight- 
ing function corresponds to a minimum intensity of 0.65 K km/s 
which means that at least 10 velocity channels show an inten- 
sity above the noise RMS. iii) Integrated intensities above some 
limit do not further increase the significance of the centroid ve- 
locities. The weighting function is unity at all points with higher 
intensities. 

Figure [10] shows the resulting A-variance spectra of the cen- 
troid map when weighting functions with different upper in- 
tensity thresholds are applied. For an easy comparison to the 
previous results we also include the graph obtained without 
any weighting. The introduction of the weighting function re- 
sults in a completely different A-variance spectrum at large 
lags. Whereas the previous computation showed a flat spec- 
trum we find now a strong increase of the spectrum above 3°. 
Unfortunately, the results only give rough guidelines for the se- 
lection of the optimum weighting function. The best continua- 
tion of the A-variance spectra from the smaller scales provided 
by the KOSMA map is obtained when the upper intensity limit 
corresponds to four times the lower limit. The two curves only 
using a narrow dynamic intensity range for the weighting func- 
tion are still heavily influenced by observational noise visible as 
increased A-variance values at small lags. On the other hand it 
is not clear whether the complete suppression of the noise ef- 
fect in the two curves with the highest upper limits is realistic, 
so that we conclude that in this example the significance of the 
centroid velocities does not further increase for integrated in- 
tensities above 2.6 K km/s. Nevertheless, a final answer to this 
question still has to come from a theoretical model for the quan- 



4. Conclusions 

In paper I we have proposed two essential improvements of the 
A-variance analysis. Here, we have tested their actual impact 
when applied to data sets characterising interstellar turbulence. 

The first improvement was the introduction of a weighting 
function for each pixel in the map. This allows us to study data 
sets with a variable data reliability across the map and to simul- 
taneously solve boundary problems even for maps with irregular 
boundaries. Maps with a variable data reliability are eventually 
obtained in most observations, either due to a local or a tempo- 
ral variability of the detector sensitivity or the atmosphere or due 
to different integration times spent for different points of a map. 
By applying the improved A-variance analysis to observed data 
we find that only the use of a significance function to weight 
the different data points allows us to distinguish the influence 
of variable noise from actual small-scale structure in the maps. 
In the analysis of intensity maps the weighting function is best 
provided by the inverse RMS in the data points. The situation 
is more complex for derived quantities, like centroid velocities, 
without a simple analytic relation between the uncertainty of the 
quantity and the observational noise. Here, in general two thresh- 
olds can be defined - a lower threshold below which all data 
have to be ignored and an upper threshold above which the sig- 
nificance of the data is not further improved by decreasing the 
noise. For the centroid velocities this means that the integrated 
line intensities between the two boundaries may serve as weight- 
ing function. 

The second improvement of the A-variance analysis is its op- 
timisation with respect to the shape of the wavelet used to fil- 
ter the observed maps. The application of different filters in the 
analysis of hydrodynamic simulations confirmed the result from 
paper I, that a Mexican-hat filter with a diameter ratio v = 1 .5 is 
well suited to resolve prominent structure scales and to measure 
the slope of the turbulent cascade, however, it turned out that the 
impact of the detailed shape of the A-variance filter is less sig- 
nificant for realistic data than for the artificial test data used in 
paper I. The turbulent structure is well resolved for a wide set of 
filters as long as one consistent filter shape is used throughout 
the full analysis of a data set. Comparing the density and ve- 
locity structure of the simulations shows a small but significant 
shift between the scale of the most prominent velocity structures, 
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created by the energy injection, and the most prominent density 
structures produced by the velocity field. 

Applying the new method to the exa mple of the dust emis- 
sion map of p Oph bv lMotte et al.l (1 19981) shows that the spatial 
scaling behaviour there can be described perfectly by a power 
law interconnecting the range of small clumps and more massive 
cores. The method can reproduce the size of the dominant cores 
and we find no indication for large-scale correlation between the 
clumps and cores in the data. The A-variance spectrum shows no 
break in the scaling behaviour between core s and condensation s 
in contrast to the mass spectrum derived by Mot te et all (Il998l) . 
The reason for the different behaviour of the two measures has 
to be topic of a future investigation. 

In the example of the analysis of the velocity structure in the 
Polaris Flar e we show that the powe r-law scaling behaviour es- 
tablished by Ossenkopf & Mac Low (2002) for the small scales 
is continued to large scales. However, a plateau in the A-variance 
spectrum around 5 pc indicates that the existing large-scale ve- 
locity gradient is not converted directly into a turbulent cascade. 
A possible explanation for this behaviour is the existence of a 
shock producing the large-scale gradient which is only converted 
into turbulent energy on the scale of the individual density fluc- 
tuations in the cloud. This scenario would be consistent with the 
affiliation of the molecu lar cloud to a large H I supershell by 
iMeverdierks et al.l (119911) . 

Combining the results from the turbulence simulation with 
the analysis of the Polaris Flare velocity structure indicates that a 
large-scale velocity field does not automatically produce density 
structures on those scales, but that a full turbulence cascade cov- 
ering density and velocity fluctuations evolves predominantly at 
seeds of primordially existing density fluctuations, which may 
have been produced by previous velocity fields on larger scales. 
When interpreting turbulent structures in interstellar clouds it 
has to be taken into account that close to the scale of the en- 
ergy injection a statistical analysis of the turbulent cascade is 
always affected by low number statistics as few density "seeds" 
may dominate the shape of the scaling relations there. A reliable 
statistics is only given on smaller scales. 
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